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The Carina arm region, containing the supernova remnant SNRG284.3-1.8, the high-energy (HE; E>100 MeV) binary 
1FGL J1018. 6-5856 and the energetic pulsar PSR J1016-5857 and its nebula, has been observed with the H.E.S.S. telescope array. 
The observational coverage of the region in very-high-energy (VHE; E>0.1TeV) 7-rays benefits from deep exposure (40 h) of the 
neighboring open cluster Westerlund 2. The observations have revealed a new extended region of VHE 7-ray emission. The new VHE 
source HESS J1018-589 shows a bright, point-like emission region positionally coincident with SNRG284.3-1.8 and 1FGL J1018.6- 
5856 and a diffuse extension towards the direction of PSR J1016-5857. A soft (r=2.7±0.5 s tat) photon index, with a differential flux 
at ITeV of N =(4.2±l.l)-10" 13 TeV" 1 cm 2 s 1 is found for the point-like source, whereas the total emission region including 
the diffuse emission region is well fit by a power-law function with spectral index r=2.9±0.4 sta t and differential flux at 1 TeV 
of N =(6.8±1.6)40~ 13 TeV^ 1 cm~ 2 s _1 . This H.E.S.S. detection motivated follow-up X-ray observations with the XMM-Newton 
satellite to investigate the origin of the VHE emission. The analysis of the XMM-Newton data resulted in the discovery of a bright, 
non-thermal point-like source (XMMU J101855. 4-58564) with a photon index of r=1.65±0.08 in the center of SNRG284.3-1.8, 
and a thermal, extended emission region coincident with its bright northern filament. The characteristics of this thermal emission 
are used to estimate the plasma density in the region as n~0.5 cm~ 3 (2.9 kpc/d) 2 . The position of XMMU J101855. 4-58564 is 
compatible with the position reported by the Fermi-LAT collaboration for the binary system 1FGL J1018. 6-5856 and the variable 
Swift XRT source identified with it. The new X-ray data are used alongside archival multi-wavelength data to investigate the 
relationship between the VHE 7-ray emission from HESS J1018-589 and the various potential counterparts in the Carina arm 
region. 
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1. Introduction 




PSR J1016-5857 



The H.E.S.S. (High Energy Stereoscopic System) collabo- 
ration has carried out observation s of the Carina arm as 
part of the Galactic Pla ne Survey (jAharonian et al.ll2006d 
l2008aHGast et alj2011f) . The observed region includes three 
potential VHE 7 -ray emitters, SNRG284.3-1.8 (MSH 10- 
53), the high spin-d own luminosity pulsar PSR J1016-5857 
(|Camilo et al.l [200lh and the Fermi Large Are a Telescope 
(Ferm i-LAT) 7 -ray binary 1FGL J1018. 6-5856 (|Abdo et al.l 
120101) . SNRG284.3-1.8 has an incomplete radio shell (Fig. 
[1]) with a non-thermal rad io spectrum and a flux density 
of (5.4±0.8) Jy at 8.4 GHz (|Milne et al.lll989ft . Evidence of 
interaction with molecular clou d s (MC ; 12 CO J=l-0) has 
been reported by iRuiz fc Mavl (fl98l using observations 
with the 1.2 m Columbia Millimeter- Wave Telescope. The 
MC content integrated over velocities in the local standard 
of rest from -14.95kms _1 to -21.45kms~ 1 traces the radio 
shell shape of the SNR. 

The shell-like SNR shows a bright, narrow filament, co - 
incident with an Ha filament (|Van der Bergh et al. I [l973h 
to the East. Optical and CO observations imply that 
SNRG284.3-1.8 is most likely a 10 4 yr old remnant of 
a type II supernova with a massive stellar progenitor 
and that it is located at a distance of ds , jy j R=2.9kpc 
(with an error of ±20%) (|Ruiz fe Mavlll986t) . The west- 
ern edge of the shell shows a so-called finger emis- 
sion region (see Fig. [T]) extended towards the direction 
of PSRJ1016-5857. This Vela-like pulsar and its asso- 
ciated X-ray pulsar wind nebula (PWN), located «35' 
away from t he geometrica l cente r of th e SNR, were dis- 
covered by ICamilo et all (|200ll . I2004D in a search for 
counterparts of the unidentified source 3EG J1013-5915 
(jHartman et al.l[l999t) with the Parkes telescope. The en- 
ergetic pulsar, formerly associated with the SNRG284.3- 
1.8, has a rotation period of 107 ms, a characteristic age 
T c =21kyr and a spin-down luminosity of 2.6-10 36 ergs^ 1 . 
Its distance can be estimated fro m the pulsar's disper - 
sion measure to be dpsR=9t?kpc dTavlor fc Cordeslll993l ) 
or dpsR«8kpc (ICordes fc Laziol 2002ft . Later observations 
with Chandra (jCamilo et al.l l2004t ) resulted in the detec- 
tion of the associated X-ray PWN, apparently located 
at the tip of the finger of the SNR radio emission. The 
PWN spectrum is well fit by a power-law model with 
photon index rx=l-32±0.25, absorbing neutral hydrogen 
column density Nh=(5.0±1.7)-10 21 cm~ 2 , and unabsorbed 
flux F .8-7kcV=2.8-10~ 13 ergcm~ 2 s~ 1 . The good match 
in both width and position angle of the PWN and the 
SNR finger suggests a possible connection between the 
two. However, the association of PSR J1016-5857 with 
SNRG284.3-1.8 is questioned by the large difference in 
their distance measurements and the large offset between 
the ce nter of the shell and the pulsar position. ICamilo et al.l 
pool argue that the uncertainties in the model for the free 
distribution of electrons in the Galaxy are large along spi- 
ral arm tangents, affecting the estimation of the distance 
inferred from the dispersion measure. A proper motion ve- 
locity of vpsR=500-(d/3kpc)(21 kyr/r c ) kms -1 has to be 
invoked to connect the pulsar with the geometric center 
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Fig. 1. The MOST-Molonglo radio map (in units of 
Jy/beam) of SNRG284.3-1.8 is displayed together with 
the position of the Fermi-LAT source 2FGL J1019.0-5856 
(white dashed line showing the 95% uncertainty in its po- 
sition) and its associated X-ray source (green dot, XMMU 
J1018554-58564). To the west of the SNR the radio finger is 
visible extending toward the direction of PSR J1016-5857 
(marked with a yellow filled star) . The PSF size of the radio 
image is 45" x 52.5". 



of the SNR. Even if the proper motion velocity is within 
the range of known pulsar velocities, this highly supersonic 
movement in the local interstellar medium (ISM) would 
create a bow shock nebula, which has not been observed. 

The Fermi-LAT collaboration has reported the detec- 
tion of a new source 1FGL J1018.6-5856 (or 2FGL J1019.0- 
5856), positionally coincident with SNR G284.3-1.8 and 
associated with t he E GRET source 3EG J1013-5915 
(lAbdo et al.l l2010l 120111) . The source reported in the 
LAT second-year source catalog shows an energy flux of 
(3.53±0.11)-10- 10 ergcm- 2 s- 1 in the 100 MeV to 100 GeV 
energy band, and a fit to a power-law spectrum yields 
a photon index of r HE =2.59±0.04. The AGILE col- 
laboration also reported HE emission from this re- 
gion and pulsed emission is marginally detected (4.8a) 
from the direction of PSR J1016-5857 (jPellizzonni et all 
l'2009h . The 7-ray emission from 1FGL J1018. 6-5856, how- 
ever, shows a periodic modulation with a period of 
16.58 ±0.02days, detec ted by Fermi-LAT ([Corbet et al.l 
l201lt lAckermann et alJl2012ft . This periodicity is taken as 
evidence that 1FGLJ1018. 6-5856 is a new 7-ray binary 
system, with an 06V((f))-type star (2MASS 10185560- 
5856459) proposed as the stellar counterpart. The spec- 
trum of the periodic source exhibits a break at sal GeV with 
best-fit values of r HE (0.1-l GcV)=2.00±0.04 and T HE (1- 
10 GeV)=3.09±0.06 and an integral energy flux above 
100 MeV of (2.8±0.1)-10" 10 crgcm" 2 s~ 1 . The binary na- 
ture of the HE source is strengthened by the detection of a 
periodic X-ray (Swift XRT) and radio (Australia Telescope 
Compact Array) compact source coincident with the 7-ray 
source. A counterpart in hard X-rays ha s also been detected 
using INTEGRAL data (|Li et al.ll20fl . 
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Fig. 2. H.E.S.S. excess image of the HESS J1018-589 region smoothed with a Gaussian of width cr=0.07° (on the left, 
a) and o=0.11° (on the right, b)). Significance contours, calculated using an ovcrsampling radius of 0.1° (in green) and 
0.22° (in white) are shown, starting at 4cr in steps of la. The black crosses (A and B) mark the best-fit position and 
statistical errors at la level for the two emission regions (see text for details). The position of 1FGL J1018. 6-5856 is 
shown with a blue dashed ellipse (95% confidence level) while the position of PSR J1016-5857 is marked with a yellow 
star. The white circles illustrate the size of the (68% containment radius) PSF smoothed with the two different Gaussian 
widths used. 



2. H.E.S.S. observations 

H. E.S.S. is an array of four VHE 7-ray imaging atmospheric 
Cherenkov telescopes (IACTs) located in the Khomas 
Highland of Namibia (23°16'18" S, 16°30'00" E). Each 
of these telescopes is equipped with a tessellated spheri- 
cal mirror of 107 m 2 area and a camera comprised of 960 
photomultiplier tubes, covering a large field-of-view (FoV) 
of 5° dia meter. The syste m works in a coincidence mode 
(see e.g. iFunk et all I2004D . requiring at least two of the 
four telescopes to trigger the detection of an extended air 
shower (EAS). This stereoscopic approach results in a high 
angular resolution of ~6' per event, good energy resolu- 
tion (15% on average) an d an efficient background rejection 
(|Aharonian et all l2006bl ). The mean (energy-dependent) 
point spread function (PSF) is estimated to be «0.1°. 
These characteristics allow H.E.S.S. to reach a sensitivity 
of «2.0-10~ 13 ph cm- 2 s _1 (equivalent to 1% of the Crab 
Nebula flux above 1 TeV) , o r better if advance d techniques 
are used for image analysis (jAcero et al.ll2009h . for a point- 
like source detected at a significance of 5a in 25 hours of 
observation at zenith. 

The Carina arm region benefits fro m deep H.E.S.S. ob- 
serva tions of the Westerlund 2 region (jAbramowski et ahl 
l2010h and an acceptance-corrected effective exposure time 
of 40 h was obtained towards the direction of SNRG284.3- 

I. 8 and PSR J1016-5857 in a multi-year observation cam- 
paign in January and March 2007, April and May 2008, 
and May to June 2009. The data set consists of scan-mode 
observations and dedicated observations in wobble-mode, 
in which the telescopes are pointed offset from the nominal 
source location to allow simultaneous background estima- 
tion (with an offset range in the total data set from 0.5° 



to 2.6°). The observations were performed in a zenith an- 
gle range from 35° to 50° and the optical response of the 
system was estimated from the Cherenkov lig ht of single 
muons as explained in lAharonian et ahl (|2006b ). 

The data have been analy zed using a multivariate 
analysis (|Becherini et al.l 12011) and c ross-checked with 
the Hillas se cond moment (lHillaslll984t) event reconstruc- 
tion scheme (|Aharonian et al.ll2006bl) and Model Analysis 
(|de Naurois fc Rollandl l2009t) . including independent cali- 
bration of pixel amplitudes and identification of problem- 
atic or dead pixels in the IACTs cameras, leading to com- 
patible results. To study the morphology of the source 
with sufficient statistics, a cut of 80 photoelectrons (p.e.) 
on the intensity of the EAS is used to reduce the data 
and produce the images, light curve and spectrum of the 
source. The energy threshold of the analysis presented here 
is E th «0.6TcV. 

Figure [U a) shows an image of HESS J1018-589 
smoothed with a Gaussian of width er=0.07°. The back- 
ground in each pixel is estimated wi th the ring back- 
ground method (jAharonian et a l. 2006b). This image shows 
a point-like 7-ray excess centered on the position of 
1FGLJ1018.6-5856 and SNRG284.3-1.8, and a diffuse 
emission region extending in the direction of PSRJ1016- 
5857. A peak pre-trials significance of 8.3<r (which corre- 
sponds to over 6 sigma post-trials) is obtained using an 
oversampling radius of 0.10° (green contours in Fig. [2]). Fig. 
[2]b) shows the same image but smoothed using a Gaussian 
of width (7=0.11° to match the significance contours calcu- 
lated with an ove rsampling radius of 0.2 2° (optimized for 
extended sources, |Ah aroman et al.ll2006cl) . which are over- 
laid in white. An extension of the VHE emission towards the 
West is observed when using the larger radius. The white 
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Fig. 3. Profile of the VHE emission along the line between 
the peak of the point-like emission and the peak of the dif- 
fuse emission, as illustrated in the inset. Fits using a single 
and a double Gaussian function are shown in dashed and 
solid lines respectively. The positions of XMMU J101855.4- 
58564 and PSR J1016-5857 are marked with dashed and 
dotted vertical lines and red and yellow stars in the inset, 
in which the significance image obtained using an oversam- 
pling radius of 0.1° is shown. 



circles in Fig. [2] illustrate the size of the PSF smoothed with 
the two different Gaussian widths used. 

To investigate the morphology of HESS J1018-589 dif- 
ferent models were fit to the un corrected excess m ap us- 
ing the Sherpa fitting package (iRefsdal et al.l [2009T) . The 
maximum- likelihood-ratio (MLR) test, performed using a 
Gaussian-shape function convolved with the H.E.S.S. PSF, 
finds a point-like source at a=10 h 18 m 59.3 s ±2.4| tat and 6=- 
58°56'10"±36" st at (J2000), marked in black in Fig. (posi- 
tion A). The sys tematic error is estimated to be 20" per axis 
(|Gillessenl I2004T) . Including a second source in the MLR- 
test results in an improvement of the test statistics by 29.8 
(or «4.4ct for a four degree of freedom % 2 -test). The best- 
fit position of the second source (position B) is located 
at a=10 h 16 m 31 s ±12| tat and <5=-58°58'48"±3' st at (J2000) 
with an extent of CT=0.15°±0.03° st at- 

The two-emission-regions hypothesis is strengthened 
when the projection of a rectangular region on the uncor- 
rected excess map along the extension of the image (with 
a width of twice the H.E.S.S. mean PSF) is evaluated. The 
fit of a single Gaussian function to the projection yields 
a x 2 /v of 21.04/13 (null hypothesis probability P=0.072) 
whereas the fit to a double Gaussian function leads to a 
better X 2 /v of 9.6/10 (P=0.476). Fig. H shows the profile 
of the rectangular region (in white in the inset). The single 
and double Gaussian function fits are shown in dashed and 
solid lines, respectively. 

The energy spectrum (shown in Fig. U) for the 
VHE excess was computed b y means of a forward- 
folding maximum likelihood fit dPiron et alJl200ll) . To es- 
timate the background the reflected background method 
was used, in which symmetric regions, not contami- 
nated b y known sources, are used to extract the back- 
ground ( Aharonian et a l. 2006b). To derive the source spec- 
trum two regions were selected, one around the point- 
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Fig. 4. VHE photon spectrum of HESS J1018-589 for a 
point-like source at position A (in blue dots and dashed 
blue line) and derived from a region of size 0.30° compris- 
ing the point-like and diffuse emission (in black dots and 
solid black line.). The residuals to the fit are shown in the 
bottom panel. 

like source and a second region of 0.30° radius ac- 
counting for the total VHE 7-ray emission centered on 
a=10 h 17 m 45.6 s and <5=-59°00'00"(J2000). The photon 
spectra in both cases arc well represented with a simple 
power-law function dN/dE=N (E/l TcV)~ r . The point- 
like source has a photon index of r=2.7±0.5 s t a t±0.2 sys 
statistically compatible with the total (0.30°) VHE 7-ray 
emission region, r=2.9±0.4 st at±0.2 sys . The normalization 
constants at 1 TeV are N = (4.2±l.l)-10- 13 TeV^cm^s- 1 
and N =(6.8±1.6)-10- 13 TeV- 1 cm^ 2 s~ 1 for the point-like 
source and the total emission region, respectively. The sys- 
tematic error on the normalizati on constant Np is estimat ed 
from simulated data to be 20% (lAharonian et al.ll2006bD . 

Finally, Fig. [5] shows the light curve of the point-like 
emission region, integra ted between .6 TeV and 10 TeV. 
The Lomb-Scargle test (|Scargle 111982ft was applied to the 
data to search for periodicity or variability without positive 
results for this data set. Fitting a constant to the integrated 
flux yields a \ 2 l v of 67.2/66, equivalent to a variable inte- 
gral flux on a run- by- run basis at the level of I.Oct. The lack 
of orbital coverage prevents any firm conclusion on variabil- 
ity of the TeV emission at the Fermi-LAT reported period 
(16.58 days). 



3. X-ray observations with XMM-Newton 

XMM-Newton observations were acquired (ID: 0604700101, 
PI. E. de Ofia Wilhelmi) to investigate the origin of the 
VHE emission region. The field was observed on the 22 nd 
August 2009, with a total integration time of 20ksec. 
The observations were centered on a=10 h 18 m 55.60 s and 
( 5=-58 o 55 , 56.8 ,/ (J20 00) and acquired with the EPIC-PN 
(IStriider et al.ll2001h and EPIC-MOS dTurner et al.ll2001ft 
cameras in full-frame mode with a medium filter in a sin- 
gle pointing. This position allowed a good coverage of the 
whole SNR structure in the EPIC cameras. 

The data were analyzed using the XMM Science 
Analysis System (SAS vll.O.oEl) and calibration files valid 
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Fig. 5. Light curve of the integral flux between 0.6 TeV and 
lOTeV in a 0.1° region centered on position A binned such 
that each point represents the data in one H.E.S.S observ- 
ing run (typically 28 min). An upper limit for 99% confi- 
dence level is shown with an arrow when the fluctuation 
in the integral flux is negative (grey points). The dashed 
horizontal line shows the mean integral flux. 

at September 2011. To exclude high background flares, 
which could potentially affect the observations, light curves 
were extracted above 10 keV for the entire FoV of the EPIC 
cameras, but no contamination was found. Therefore the 
full data set was used for the image and spectral analy- 
sis. To create images, spectra, and light curves, events with 
FLAG=0, and PATTERN=12 (MOS) and 4 (PN) were se- 
lected. Hereafter clean event files in the 0.3 to 8keV energy 
band are used. 

Images combining the different EPIC instruments (see 
Fig. [6] a) and b)), vignetting-corrected and subtracted for 
particle induced and soft proton background, were pro- 
duced using the ESAS analysis package (integrated in SAS). 
In the 0.3 keV to 2keV energy range a bright source is de- 
tected (see Fig. [5] a)). The radial profile of this central 
source was derived from the three EPIC cameras and fit 



with the corresponding PSF, confirming its point-like na- 
ture within the instrument angular resolution and observa- 
tion sensitivity. This source is surrounded by diffuse emis- 
sion extending up to the radio shell. A strong enhancement 
of the diffuse emission is visible just downstream of the 
radio and Ha filament (see Fig. [5] a) and d)). At higher 
energies (2 to 8keV, see Fig. [5] b)), the diffuse emission 
is strongly suppressed suggesting a thermal nature of the 
emission pervading the SNR, the only significant feature 
being the bright point-like central source. 

This source, dubbed XMMU J101855.4-58564, is lo- 
cated at a=10 h 18 m 55.40 s and c5=-58°56'45.6" (J2000) with 
statistical error of ±0.25" in each coordinate (derived us- 
ing the SAS task e detect). The p o sition is compatible with 
the one derived by iPavlov et al.l (|2011| ) The photon spec- 
trum of the point-like source was derived integrating over 
a 20" circle around the fit position and the background 
was estimated from a circle of 40" located in the vicinity. 
The spectrum is well fit (x 2 ed =0.97, z^=159) by an absorbed 
power-law function in the 0.5 to 7.5 keV energy range, with 
a photon index of T x — 1.66±0.11 s t a t and an integrated flux 
of F2_iokcV=(6.5±0.7 s tat)T0" 13 ergcm- 2 s- 1 . The absorp- 
tion column density Nh is (6.6±0.8 sta t)T0 21 cm~ 2 , sup- 
porting a Galactic origin of the source. Other models such 
as a black-body model give fits that are statistically in- 
adequate. Fig. [7] shows the measured spectrum for MOS1, 
MOS2 and PN (in black, red and green respectively). 

Archival 2MASS (Two Micron All Sky Survey) data 
of the region show a bright star (with magnitudes 
J=10.44±0.02, H=10.14±0.02 and K=10.02±0.02) dubbed 
2MASS 10185560-5856459, located at a=10 h 18 m 55.6 s and 
<5=-58°56'46" (J2000), 1.3" away from XMMU J101855.4- 
58564, the likely counterpart in the binary system. This 
source also appears in the USNC0 catalog with magnitudes 
B=12.76 and R=11.16. A distance of d*=5.4±^kpc to 
th e 2MASS star has b een estimated through photometry 
bv lNapoli etHl (|201lh . The position of XMMU J101855.4- 
58564 is also in agreement with the variable com- 
pact object detected by Swift XRT at a=10 h 18 m 55.54 s 
and <$=-58°56'45.9" (J2000) and the Fermi-LAT source 
1FGL J1018. 6-5856 dCorbet et al.l[20Tl . 

To the North-East of XMMU J101855.4-58564 a faint 
extended emission region located just downstream of the 
radio and Ha filament of the remnant is visible at low en- 
ergy (E<2kcV) (Fig. E] a) and d)). To extract the X-ray 
spectrum, the background was m odeled using the ESA S 
software following the approach of ISnowden et all (|2004f ) . 
This background model is subsequently used to fit the sig- 
nal region. The final shell spectrum is well represented by 
an absorbed non-equilibrium ionization (PSHOCK) ther- 
mal model with a temperature of kTps0.5 keV and a column 
density of 8T0 21 cm~ 2 (in Fig. [gj. The normalization factor 
A, defined as 



A(cm- 5 ) = 10" 



n 2 dV/47rd 2 



(1) 



is 1.5T0 -3 cm -5 . From this value, assuming a spherical vol- 
ume (V) corresponding to the 2' source extraction region 
(see Fig. [Sk)) and a fully ionized gas (n e =1.2-n), a plasma 
density can be derived to be n«0.5cm _3 (2.9kpc/d) 2 . 



http:/ /www.usno.navi/mil/USNO/ 
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Fig. 6. Multi- wavelength view of the HESS J1018-589 region. On the top left a) XMM-JVewfon exposure-corrected counts 
map centered on SNRG284.3-1.8 in the 0.3 to 2keV energy band in color scale combining the three EPIC cameras. The 
image (in excess per arcmin 2 ) is in logarithmic scale and it was adaptively smoothed to emphasize the extended thermal 
emission. The white-dashed contours show the radio image of SNR G284.3-1.8 from MOST-Molonglo observations at 
843 MHz. The extraction region to evaluate the diffuse X-ray emission is shown in green. On the top right b) XMM- 
Newton exposure-corrected counts map with the EPIC cameras at high energies in the 2 to 8 keV energy band. The white 
cross and yellow circle mark the best fit position of the H.E.S.S. point-like emission and the 95% confidence level of the 
Fcrmi-LAT source 2FGL J1019. 0-5856. On the left bottom c) The significance contours for the two different ovcrsampling 
radii for HESS J1018-589 arc shown in green (0.1°) and white (0.22°). The best-fit positions (white cross) of the VHE 
7-ray, superimposed on the MOST-Molonglo radio map are displayed, together with the position of PSR J1016-5857 
(yellow star) and 2FGL J1019. 0-5856 (in dashed- yellow) . On the right bottom d) A zoom of the radio shell (in white) 
is shown together with the contours of the thermal emission in green and Ha observations from the SupcrCOSMOS 
H-alpha Survey (in intensity units). 



4. Discussion 

A new VHE 7-ray source, HESS J1018-589, has been 
discovered in the vicinity of the shell-like remnant 



SNRG284.3-1.8. The extended TeV emission region coin- 
cides with SNRG284.3-1.8 and PSR J1016-5857, both vi- 
able candidates to explain the observed VHE 7-ray emis- 
sion. This emission region is also positionally compatible 
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Fig. 7. XMM-iVewtorc/PN/MOS spectrum of the point- 
like source XMMU J101855. 4-58564, fit with an absorbed 
power-law model in the 0.5 to 7.5 keV energy range. The 
bottom panel shows the residuals to the fitted spectra. 
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The new HE 7-ray binary 1FGL J1018. 6-5856 shares 
many characteristics wit h the VHE 7-ray binary LS 5039 
()Aharonian et al.l|2005al) . The XMM- Newton observations 
presented here reveal a bright non-thermal point-like 
source, XMMU J101855.4-58564, in the center of the SNR 
and compatible with the position of the binary system and 
the H.E.S.S. point-like emission. The X-ray photon spec- 
trum resembles that of pulsars, with a photon spectral in- 
dex of 1.67 and a column density of (6.6±0.8)-10 21 cm -2 , 
compatible with that from the thermal emission region co- 
incident with the bright radio filament (7.9T0 21 cm -2 ). The 
position of XMMU J101855. 4-58564 centered on the SNR 
and the similar distance (also compatible with the dis- 
tance to the associated 2MASS star) suggest a physical 
association between the two objects (rather than an as- 
sociation between the SNR and PSR J1016-5857), namely 
that the compact object within the binary system was the 
stellar progenitor for the type II SN explosion (see e. g. 
iBosch-Ramon fe Barkovl 1 2 llh . No extended emission or 
putative PWN around the point-like source is observed in 
the present data at the level of the XMM-Newton observa- 
tion sensitivity. The light curve of the X-ray emission docs 
not show any indication of variability or periodicity on short 
time scales, for the time resolution of the MOS and PN 
cameras in full-frame mode (2.6 s and 73.4ms respectively) . 
XMM U J101855.4-58564 has been associated (|Pavlov et all 
[20T1 with the Swift XRT source, which itself has been 
identified as a counterpart to the HE binary system. The 
best-fit position of HESS J1018-589 is compatible within 
less than la with both the position of XMMU J101855.4- 
58564 and the variable HE source, whose position has been 
determined accurately using timing analysis. The spectral 
type of the possible companion star, 06V((f)) is similar to 
the one in the VHE 7-ray binary LS 5039. 

In a binary scenario composed of a massive star and 
a pulsar or a black hole, modulated VHE 7-rays can 
be produced by different mechanisms, namely inverse 
Compton emission or pions produced by high- energy pro- 
tons interact i ng with the stella r wind dBednarekl [200(: 
Dubusl 120061: iKhaneulvan et all 120081: iKirk et all ll99S 



Fig. 8. XMM- Newton /PN /MOS spectrum of the extended 
emission region (in green in Fig. El a)) fit with an absorbed 
non-equilibrium ionization (PSHOCK) thermal model. The 
bottom panel shows the residuals to the fitted spectra. 



with the new HE 7-ray binary 1FGL J1018. 6- 5856 reported 
recen tly by the Fermi-LAT collaboration (jCorbet et al.l 
120111) . Different possible scenarios in the context of multi- 
wavelength observations, including new results derived 
from XMM-Newton observations, arc discussed in the fol- 
lowing. 

Although the spectral characteristics and light curve of 
the H.E.S.S. source do not yet allow a firm identification of 
the origin of the VHE 7-ray emission, the morphology of 
the source is considered to clarify the situation. Two dis- 
tinct emission regions are detected in the H.E.S.S data, 
one point-like emission region (A) located in the center 
of SNRG284.3-1.8 with a centroid compatible with the 
95% confidence contour of 1FGL J1018. 6-5856, and a dif- 
fuse emission region (B) extending towards the direction 
of PSR J1016-5857, with its centroid compatible with the 
position of the pulsar. 



Sierpowska-Bartosik fe Torres! l2008h . Similar to previously 



detected VHE binar y systems such as LS 5 039, LSI +61 303 
(I Albert et al.l I2009D and PSRB1259-63 (jAharonian et al.l 
l2005bl ) periodic emission is found at HE. Despite the simi- 
larities with other VHE binaries and in particular LS 5039, 
and the good positional agreement with the Fermi-LAT 
source, the association of the H.E.S.S. source (A) with 
1FGL J1018. 6-5856 is still uncertain. No flux variability has 
been observed yet and the Lomb-Scargle test does not re- 
cover the Fermi-LAT reported 16.58 days modulation of the 
HE signal. Nevertheless it should be noted that the non- 
detection might be due to possible contamination from the 
neighboring diffuse emission, statistics of the data set, and 
and inadequate time sampling of the orbit. 

HESS J1018-589 (A) is also coincident with the 
SNRG284.3-1.8. SNRs are believed to be sites of parti- 
cle acceleration up to at least a few tens of TcV. Two types 
of VHE 7-ray emission associated with SNRs have been 
discovered with IACTs, VHE 7-r ay emission from shell- like 
SNRs suc h as RXJ1713. 7-3 947 (jAharonian et al.ll2007l) or 
SN 1006 (|Acero et al.ll2010l) . in which in general the VHE 
morphology is in good agreement with the synchrotron X- 
ray emission; and VHE 7-ray radiation which seems to 
originate through proton-proton (p-p) interaction of cos- 
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mic rays (CR) accelerated in the SNR interacting with 
local MCs in the vic i nity, s uch as W28 (jAharonian et al.l 
l2008bf) . iRuiz fe Mavl (|1986| ) reported evidence of interac- 
tion of SNRG284.3-1.8 with embedded MC, constraining 
the distance to the SNR to «2.9kpc. The analysis per- 
formed in the context of this work of public 12 CO (J =1- 
0) data from the CfA 1.2 m Millimeter- Wave Telescope 



(|Dame et al.l 20011) yields an est imate of the MC mass of 



3-10 d M 



©• 



Ruiz fc Mavl (|1986l) also reported on optical 
observations in the direction of SNRG284.3-1.8 and asso- 
ciated a bright optical filament (Fig. E]d)) coincident with 
the brightest shell structure to the North-East, indicating 
collisional excitation of the ISM, swept by the expanding 
SNR shock wave. The observed MCs, if indeed physically 
associated with the SNR, could provide enough target ma- 
terial to explain the VHE emission in a scenario in which 
the 7-ray are produced via p-p interaction. However, con- 
trary to some other SNRs at VHE, the emission detected 
with H.E.S.S. does not match the shell-type morphology 
within the present statistics. 

The diffuse emission detected with H.E.S.S. (B) extends 
towards the direction of PSR J1016-5857. PSRJ1016- 
5857 was detected as a bright EGRET source and 
pulsed e mission at HE has been reported by the 
AGILE ijPellizzonni et al.| 
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and the Fcrmi-LAT 
(|Saz Parkinson et al. 2010 ) collaborations. With a spin- 
down luminosity of 2.6T0 36 ergs -1 , the radio, HE and X- 
ray pulsar PSR J1016-5857 is energetic enough to power 
the entire H.E.S.S. source, assuming a dispersion-measured 
estimated distance of 9 kpc. In this scenario, particles are 
accelerated in the wind termination shock and produce 
VHE 7-ray emission by inverse Compton (IC) processes 
as they propagate away from the pulsar. As a result of 
the interactions of relativistic leptons with the local mag- 
netic field and low-energy radiation, non-thermal radia- 
tion is produced up to »100Te V (for a recent review see 
Ide Jager fc Diannati-Atal |2008| ). Assuming a distance of 
9 kpc, the total VHE luminosity in the 1 to 10 TeV energy 
range is 9.7T0 33 -(d/9 kpc) 2 ergs" 1 , implying a maximum 
conversion from rotational energy into non-thermal emis- 
sion with efficiency 0.4%, with similar feat ures to other well- 
establi shed VHE PWNe, su ch as Vela X (|Aharonian et al.l 
I2006dh or HESS J1026-582 (lAbramowski et al.ll2010D . The 
associated X-ray nebula has been detected with Chandra in 
the 0.8 to 7keV energy range with a size of 3'. The differ- 
ent size of the VHE and X-ray nebula, «20-(d/9kpc) pc and 
«8-(d/9kpc) pc respectively, can be easily accommodated 
in a relic nebula scenario, and explained by the different 
energies (and hence cooling times) of the electron popu- 
lation emitting X -rays and VHE 7-rays i 
HESSJ1825-137 (lAharonian et al.lhoOBal) 
netic field of the order of a few /iG. 



as seen, e. g. m 
for a low mag- 



5. Conclusions 

A new VHE 7-ray source dubbed HESS J1018-589 has been 
detected with the H.E.S.S. telescope array with a signif- 
icance of 8.3cr. The complex VHE morphology and faint 
VHE emission prevent a unequivocal identification of the 
source given the presence of several possible counterparts. 
The H.E.S.S. source seems to be composed of two emis- 
sion regions but the statistics are still too low to make firm 
conclusions about the origin of those. 



Several counterparts are discussed using energetics ar- 
guments as to the possible origin of either part or all of the 
emission. In a SNR/MC scenario, SNRG284.3-1.8 could 
partially explain the VHE 7-ray emission via p-p interac- 
tions with the associated MCs. However the fact that the 
VHE emission does not trace either the irregular shell or 
the cloud morphology disfavors SNR G284.3-1.8 as the only 
counterpart. 

The morphology and good positional agreement be- 
tween the H.E.S.S. best-fit position A and the new Fermi- 
LAT binary 1FGL J1018. 6-5856 suggest a common origin. 
The analysis of the XMM- Newton observations revealed 
a non-thermal point-like source, XMMU J101855.4-58564, 
with photon spectral index of 1.67, similar to the compact 
object found in LS 5039. Likewise, the spectral class of the 
massive star companion listed in the 2MASS and USNO 
catalogs is similar to the one in LS5039. However, no vari- 
ability has been found in the H.E.S.S. light curve. A dedi- 
cated observation campaign at VHE should help to clarify 
whether or not the two sources are indeed associated. 

The energetic pulsar PSR J1016-5857, also recently de- 
tected in Fermi-LAT and AGILE data, and its X-ray neb- 
ula seem the most likely candidate to power the extended 
VHE 7-ray source, given the high spin-down luminosity 
(2.6T0 36 ergs -1 ) and X-ray nebula, which implies a pop- 
ulation of high energy electrons able to up-scatter soft pho- 
ton fields to VHE. The estimated age of the pulsar (21 kyr) 
would also explain the large size o f the VHE nebula, simila r 
to other systems such as Vela X (jAharonian et al.ll2006dl) . 

Finally, XMM- Newton observations also revealed ther- 
mal emission behind the brightest synchrotron part of the 
radio shell of SNR G284.3-1.8, which might be associated 
with shock heated interstellar matter. The column density 
is statistically compatible with the one derived from the 
direction of XMMU J101855.4-58564. The similar column 
density and the position of the pulsar candidate with re- 
spect to the center of the SNR could indicate a common 
origin, where XMMU J101855. 4-58564 is interpreted as the 
pulsar left behind after the supernova explosion. 
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